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Abstract: The commercial use of transgenic cotton expressing an insecticidal protein gene from Bacillus thuringiensis (Bt) began in 1998 in 
northern China. Resistance management is a major concern for the sustainable use of Bt cotton. With our understanding of the cropping and 
ecological system in northern China; we developed a simulation model to forecast adaptation of the cotton bollworm, Helicoverpa armigera 
( Hübner), to Bt cotton. The model incorporated population genetics, biological and operational factors in a multiple cropping system. Vari- 
ables of reproductive fitness on different host plants; dominance of resistance alleles, expression level of Bt toxin in cotton and use of insecti- 
cides in Bt cotton fields were included in the model. In the typical cropping system of northern China, the expected life of Bt cotton are seven 
years if all cotton is Bt cotton, and ten years if only spring planted cotton (about 70% total cotton area) is Bt cotton in northem China besed 
on the model. The life expectancy decreases quickly with increases in initial frequency of resistance allele. dominance of the resistance gene» 
and the percent area of Bt cotton planted. The results also showed that supplemental control is essential on Bt cotton when the expression of Bt 


toxin declines. Because of the lack of a consistent high dose and the complexity of managing individual fields of Bt cotton by many different 


growers in the same region, the use of the high-dose/refuge strategy in northern China will be difficult to achieve. 
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The cotton bollworm, Helicoverpa armigera 
(Hübner); is the most important cotton pest in China. 
Long-term application of pyrethroids resulted in serious 
resistance problems and reduced control efficacy (Fan et 
al.» 1997). Transgenic cotton» developed in the US and 
China and expressing an insecticidal protein from Bacillus 
thuringiensis ( Bt), provided a powerful new tool after 
commercial use began in northern China in 1998. The 
acreage of Bt cotton grown in China was about 80 000 ha 
in 1998 and over 0.3 million ha in 1999, and greater 
than 90% of the cotton fields in Shandong and Hebei 
provinces in northern China cotton belts were Bt cotton in 
2000 (Zhao et al.» 2000). A survey of 282 Bt cotton 
farmers in northern China demonstrated significant in- 
creases of production efficiency and decreased health risks 
to farmers (Huang et al.» 2002). To preserve these ben- 
efits a resistance management program must be implement- 
ed so that insects do not rapidly become resistant to the Bt 


cotton since some insects have already demonstrated their 


ability to develop resistance to Bt toxins Shelton et al.» 
2002>. 

The high-dose/refuge is the only resistance manage- 
ment strategy currently available for Bt transgenic crops» 
and the US growers must chose one of the structural refuge 
options as required by the US EPA CShelton et al.» 
2002). To use Bt plants expressing a high-dose of the 
toxin in combination with a refuge would conserve suscep- 
tible alleles within the population. The high dose is 
required through the entire crop’ s growth to ensure that 
homozygous and heterozygous susceptible insects are con- 
trolled and the refuge strategy must be considerate of the 
ability of insects in the refuge to mate with individuals 
carrying resistant alleles» thereby decreasing the proba- 
bility of resistant phenotypes. Bt cotton provides excellent 
control of Heliothis virescens in the US (Gould and Ta- 
bashnik» 1998), but Helicoverpa zea in the US and H. 
armigera in China are less sensitive to the available Bt 


toxins and showed significant survival on Bt cotton in the 
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fields (Stone and Sims» 1993; Zhao et al.» 1998, 
2000). The survival of H . zea on Bt cotton in the US in- 
dicated that the “high-dose” approach is not valid for this 
species (Gould and Tabashnik, 1998). 

Pest resistance management (PRM? strategies for 
transgenic crops have been examined in several models 
(Caprio, 1994; Tabashniky 1994; Onstad and Gould: 
1998). The fundamental goal is to minimize exposure of 
pests to the toxin in space and time. The models reported 
that expression level of Bt toxin and the insect’ s sensitiv- 
ity are key factors in resistance development. Additional- 
ly» it is important to have sufficient susceptible alleles to 
dilute the frequency of resistance alleles» and that there 
be random mating of insects carrying each type of allele. 
Models and experimental data also showed that separate 
but adjacent refuges might be superior to other strategies 
for insects which can move between plants in their larval 
stage (Shelton et al.» 2002). We developed a simulation 
model that includes reproductive fitness of H. armigera 
on different host plants» dominance of resistance alleles. 
expression level of Bt toxin in cotton» and use of insecti- 
cides in Bt cotton fields» under ecological situations in the 


northern China cotton belt. 
1 Materials and Methods 


1.1 Model design and realization 

Development of our model are based or the following 
key assumptions. (1) Resistance of H. armigera to Bt 
toxin is controlled by a single autosomal locus ©G) with 
two alleles (Gp and Gs)» and the degree of dominance is 
h. Thus, there are three genotypes (G,Gs> susceptible: 
GaGs, heterozygous: and GaGa: resistant). The initial 
resistance allele frequency is p. (2) There are n types of 
host-plant fields planted. Wheat» cotton and corn are ma- 
jor host plants of H. armigera. Movement of larvae be- 
tween different kinds of fields is negligible. and genera- 
tions do not overlap. Mating occurs randomly among 
adults of different genotypes and different fields. Egg dis- 
tribution depends on the preference of H. armigera 
adults on different host plants. (3) Population survival is 
assumed to be density-independent. Genotype survival 
rates are defined as L for GpGp-> K for GsGs> and AL + 
(1- h?) K for GRG, CAlstad and Andow, 1995). The 
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reproductive disadvantage of resistant genotypes on non- Bt 
host plants is negligible. (4) Survivors feeding on differ- 
ent types of non-Bt crops have the same fertility rate. (5) 
Proportion of each genotype in each generation can be cal- 
culated according to the Hardy-Weinberg equation. 
1.2 Model developing progress from N generation to 
N+1 generation 
The three genotypes are designated as G, CGG)» 
G (GaGs) and Gs ©GpGa). F; Ci = 1~ n) denotes 
the ith type of field. Assuming total crop area is denoted 
by A> the percentage of F, is denoted by A;. T, denotes 
the adults oviposition probability upon the plants in F; . 
D denotes average population density of all fields in the 
egg stage. D, denotes population density on F;. The fre- 
quencies of Gp and Gs are given by p and q where q = 
1- p. The frequency of G; Cj =1~3D is denoted by P, 
(j =1~ 3). The population density of genotype G; in F; 
is denoted by D; . 
Dy = P; DACA,T,/ >) CA,T,D/CAA, ) 
= DPTJ/S)CAT (i=1~n,j=1~ 3) 
E; ©i=1~ n> j =1~3) denotes the natural survivor- 
ship of the genotype G; in F, as influenced by environ- 
mental factors. Sy Ci=1~ n» j = 1~ 3) denotes the 
relative survivorship of the genotype G, in F; as influenced 
by a supplemental insecticide and other operational factors 
Cineluding light trap to kill moths). W; Ci =1~ n) de- 
notes the effective fitness of the genotype G, in F, . 
W, = Ey Ge Ln 
Ni (i=1~n, j=1~3) denotes the number of surviv- 
ing adults of the genotype G; in F; . 
Ny = CAAD; )W; = DAP,A;T.W, 1 >) CAT, 
Giz lwn jf=21~3) 
p' denotes the frequency of the resistant allele in the egg 
stage of the next generation and is independent of D and 
A. 


p = > Na + No2XJ > Ny) 
= >) CPsAT Ws + P,A,T;Wo!2)/ 
>) (PAT 
CGi=1l~n j=1-~3) 
1.3 Initializing parameters 


In northern Chinas H. armigera can complete four 
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generations per year. The larvae of the first generation 
feed mainly on wheat. The larvae of the second generation 
feed on several crops planted in the spring» but the domi- 
nant host plant is cotton. The third and fourth generations 
feed mainly on crops planted in the summer, including 
cotton» com» peanut» soybean» and other vegetables CShi 
et al.» 1995). The percentage area of each crop is set 
according to survey data from Handan City, Hebei Prov- 
ince. We assume that Bt cotton is planted in the spring 
Cabout 70% of total cotton area) and non-Bt cotton is 
planted in the summer. Other crops Cnon-Bt) are consid- 
ered as refuges of H. armigera because they are planted 
adjacent to cotton and most of the fields are very small 
(< 1 ha). The reproductive fitness of H. armigera on 
wheat and cotton planted in the spring was calculated ac- 
cording to the results of Dai et al. (1991). The repro- 
ductive fitness on com, cotton planted in the summer and 
other crops was calculated according to the results of Shi 


et al. (1995). The oviposition performance on each crop 
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was calculated according to the results of Shi et al. 
(1995). 

Because we do not have a H. armigera strain highly 
resistant to Bt» we assume that the characterization of Bt 
resistance in H. armigera was similar to that in H. 
virescens (Gould et al.» 1995). The dominance of resis- 
tance allele A was about 0.3. The initial frequency of 
resistance allele is assumed to be 0.001. Expression of Bt 
toxin in cotton begins to decline after flowering and some 
H. 
tissuess (Zhao et al., 1998). The relative fitness of 


armigera larvae could survive on flowers and other 


Gs and Gn Gs are different in each generation. The de- 
tails are listed on Table 1. The relative fitness of Gn Gm 
was 1.0 in any situation. Assuming larval mortality 
caused by chemical insecticides is 80% on non-Bt cot- 
ton» the relative fitness of each genotype was 0.2. In oth- 
er non-Bt crop fields» each genotype has a relative fitness 


of 1.0 because H. armigera was unmanaged. 


Table 1 Initializing parameters of the typical cotton cropping system in northern China 














Field First generation Second generation Third generation Fourth generation 
Genotype 
type Aj T Egy Ss 4 T; Ey Si Ai T; Ej Si Ai T, Ey Si 
Fi Gy 0 一 一 一 0.06 0.84 0.00059 0 0.06 0.2 0.0031 0.01 0.06 0.22 0.0018 0.01 
G = 一 0.00059 0.0001 0.0031 0.31 0.0018 = 0.31 
Gs 一 一 0.00059 1.0 0.0031 1.0 0.0018 1.0 
Fy G 0 一 一 一 0 0.84 0.00059 0.2 0 0.2 0.0031 0.2 0 0.22 0.0018 0.2 
G = = 0.00059 0.2 0.0031 0.2 0.0018 0.2 
Gs = 三 0.00059 0.2 0.0031 0.2 0.0018 0.2 
Fy Gi 0 - = = o 0 0.34 0.0025 0.002 0 0.4 0.0016 0.01 
G - - = - 0.0025 0.3 0.0016 0.31 
Gs 一 一 一 一 0.0025 1.0 9. 916 1.0 
Fa G 0 一 一 一 0 0.03 0.34 0.0025 0.2 0.03 0.4 0.0016 0.2 
G2 一 一 一 一 0.0025 0.2 0.0016 0.2 
G3 一 Š - 一 0.0025 0.2 0.0026 0.2 
Fs G 0 一 一 一 0 一 一 一 0.26 0.31 0.00055 1.0 0.26 0.2 0.00041 1.0 
G2 一 三 一 一 0.00055 1.0 0.00041 1.0 
Gs 一 三 一 一 0.00055 1.0 0.00041 1.0 
Fe Gy 0.41 1.0 0.0007 1.0 0.02 0.16 0.00021 1.0 0.22 0.15 0.0005 1.0 0.22 0.144 0.00061 1.0 
G 0.0007 1.0 0.00021 1.0 0.0005 1.0 0.00061 1.0 
Gs 0.0007 1.0 0.00021 1.90 0.0095 1.0 0.00061 1.0 
Fj: Bt cotton planted in spring; Fo: non-Bt cotton planted in spring; F3: Bt cotton planted in summer; Fy: non-Bt cotton planted in summer; Fs: corn planted in summer; 


Fs: other field 
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1.4 Model programming and sensitivity analysis 
The model is implemented as a C program for DOS. 
We assumed that Bt cotton would lose its efficacy when 
the frequency of resistance allele increased to 50% . The 
expected life of Bt cotton was defined as the number of 
years required for the frequency of resistance allele to in- 
crease to 50% from the initial value. Sensitivity analysis 
was made by running the program when the initial value of 
each of certain critical variables C initial resistance fre- 
quency» gene dominance, Bt cotton area» and chemical 


control) changed. 
2 Results 


2.1 Expected commercial life of Bt cotton and resis- 
tance development of H. armigera 


In typical conditions as in Table 1; the model pre- 


0.8 
U.6 


OA 


Frequency of resistance 
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dicted that the frequency of resistant allele will rise to 0.5 
after 38 generations» i. e., approximately 10 years of 
widespread use of Bt cotton Csince 2000) in northern Chi- 
na (Fig. 1). 
2.2 Sensitivity analysis of the model 

Results of sensitivity analysis indicated that the ex- 
pected life of Bt cotton will decrease quickly with the in- 
crease of initial frequency of resistance allele (Fig. 2A), 
gene dominance (Fig. 2B), and area percentage of Bt 
cotton (Fig. 3). If all cotton fields are Bt cotton» the ex- 
pected life may be only 7 years (26 generations) (Fig. 
3). Chemical insecticide applications in Bt cotton will 
prolong the life of Bt cotton and delay resistance develop- 
ment of the pest to Bt cotton (Fig. 4) because all geno- 
types» including homozygous resistant genotypes and het- 


erozygotes, will be killed. 





Fig. 1 
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Fig. 2 Sensitivity analysis on initial frequency of resistance allele (A) and gene dominance (B) in H. armigera 
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2.3 Interaction between dominance of resistance al- 
lele and expressing level of Bt toxin 

The expression level of Bt toxin in Bt cotton is pre- 
sented as the relative fitness of heterozygotes on Bt cotton 
(Fig. 5). The arrow indicates the point at which all sus- 
ceptible individuals are killed on Bt cotton when the rela- 
tive fitness of heterozygotes was equal to h (Fig. 5). In 
general, there were many more susceptible individuals 
than resistant individuals or heterozygotes in the early 


stage. At a low-level expression of Bt toxin which is not 
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Fig. 3 Sensitivity analysis on the 
regional proportion of Bt cotton 
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very harmful to any genotypes» a Bt cotton field with sus- 
ceptible survivors could be considered an “inside refuge”. 
At a moderate expression level that kills all susceptible 
individuals» the expected life of Bt cotton will be the 
shortest» with loss of an “ inside refuge” and large num- 
bers of surviving resistant individuals and heterozygotes. 
At a high expression level that could kill all heterozygotes 
and susceptible individuals, the expected life of Bt cotton 


will increase with the increase in of heterozygotes mortali- 


ty. 
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Fig. 4 Sensitivity analysis on the 


chemical control in Bt cotton field 
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Fig. 5 Sensitivity analysis of expressing level of Bt toxin in cotton 
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3 Discussion 


Due to the declining level of expression of the Bt 
toxins when the third and fourth generation of H. armig- 
era occur they will not be 100% controlled by Bt cotton 
(Zhao et al.» 1998). The lack of a “high dose” in cur- 
rent Bt cotton cultivars for H. armigera and the small 
scale production systems of cotton in northern China indi- 
cates that the “high dose/refuge” resistance management 
strategy is probably not feasible for Bt cotton in most cot- 
ton belts in China (Zhao et al.» 2000). There is also no 
mandatory requirement for a structured refuge by the regu- 
latory authority in China because of complexity of manag- 
ing individual small fields of Bt cotton by many different 
growers in the same region. The expected life of Bt cotton 
based on our modeling study is about 7 ~ 10 years if 70% 
~ 100% cotton planted is Bt cotton, indicating that alter- 
native resistance management strategies are needed for use 
by millions of cotton growers in northern China. 

Theoretical models suggest that pyramiding two dis- 
similar toxin genes in the same plant has the potential to 
delay the onset of resistance much more effectively than 
single toxin plants released spatially or temporally 
(Roush, 1998). We have shown that transgenic tobacco 
expressing two insecticidal toxin (Bt and cowpea trypsin 
inhibitor’ CpTI) genes could delay resistance develop- 
ment by H. armigera compared with Bt tobacco (Zhao et 
al.» 1999). Transgenic cotton expressing both Bt and 
CpTI genes developed in China (Guo et al.» 1999) went 
into commercial use in 2000 and the efficacy on H. ar- 
migera was significantly higher than one Bt gene cotton 
(Zhao et al.» 2000). Stewart et al. ©2001) demonstrat- 
ed that Bollgard II cotton that produce two Bt gene prod- 
ucts CCryl Ac and Cry2Ab> has greater insecticidal activi- 
ty than the original Bt cotton expressing only the CrylAc 
toxin. 

Supplemental chemical insecticides could be used to 
control more resistant genotypes and heterozygotes in Bt 
cotton fields when the Bt toxin expression declines. Our 
modeling results proved that the expected life of Bt cotton 
can be doubled if 90% H. armigera larvae on Bt cotton 
were killed by insecticides. Different management tactics 


are needed for different cotton regions in China because of 


学 报 45 卷 第 2 期 


greatly different ecological situations among them. Model 
modification and validation for this study are needed based 


on new data of Bt resistance in H. armigera . 
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tobacco expressing two insecticidal genes to delay resistance development 


华北 地 区 棉铃 虫 对 转 Bt 基因 抗 虫 棉 
抗 性 适应 的 模拟 模型 


ee, REA, HAM 


《中 国 农业 科学 院 植物 保护 研究 所 , JERR 100094) 


摘要 : 通过 对 华北 地 区 耕作 制度 和 生态 系统 的 了 解 ， 在 充分 考虑 种 群 遗 传 学 、 生 物 学 和 人 为 操纵 因子 等 三 大 因素 的 基础 
上 ， 建 立 了 一 个 预测 棉铃 虫 对 转 Ht 基因 抗 虫 棉 抗 性 适应 的 模拟 模型 。 在 华北 地 区 典型 的 耕作 制度 下 ， 如 果 所 有 棉田 均 为 Bt 
棉 ， 则 Bt 棉 的 预期 寿命 为 7 年 ， 如 果 只 有 春播 棉 为 Bt 棉 《〈 约 占 棉田 总 面积 的 70% )， 则 其 寿命 为 10 年 。 模 型 的 灵敏 度 分 析 
表明 ,Bt 棉 的 使 用 寿命 随 抗 性 基因 的 显 性 度 、 初 始 抗 性 频率 、Bt 棉 所 占 比 例 等 因素 的 增长 而 迅速 缩短 。 当 Bt 棉 表 达 的 杀 虫 
蛋白 量 恰好 全 部 杀 死 敏感 基因 型 (GsGs) 个 体 时 ，Bt 棉 的 预期 寿命 最 短 。 由 于 国外 采用 的 “高 剂量 /庇护 所 ” 抗 性 治理 策略 
不 适用 于 棉铃 虫 及 华北 棉 区 的 耕作 制度 ， 我 国 需要 加 强 对 其 它 抗 性 治理 措施 《如 转 双 基 因 抗 虫 棉 ) 的 研究 与 应 用 。 
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